p-wave superfluid and phase separation in atomic Bose-Fermi mixture 
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We consider a system of repulsively interacting Bose-Fermi mixtures of spin polarized uniform 
atomic gases at zero temperature. We examine possible realization of p-wave superfluidity of 
fermions due to an effective attractive interaction via density fluctuations of Bose-Einstein con- 
densate within mean-field approximation. We find the ground state of the system by direct energy 
comparison of p-wave superfluid and phase-separated states, and suggest an occurrence of the p- 
wave superfluid for a strong boson-fermion interaction regime. We study some signatures in the 
p-wave superfluid phase, such as anisotropic energy gap and quasi-particle energy in the axial state, 
that have not been observed in spin unpolarized superfluid of atomic fermions. We also show that a 
Cooper pair is a tightly bound state like a diatomic molecule in the strong boson-fermion coupling 
regime and suggest an observable indication of the p-wave superfluid in the real experiment. 

PACS numbers; 03.75.Ss,03.75.Mn 



I. INTRODUCTION 

Magnetically tunable Feshbach resonances have 
opened up a new field of research in the physics of ul- 
tracold atomic gases that exhibits exciting phenomena. 
In two-component fermi gases, for instance, the conden- 
sate of atom pairs in the BCS-BEC crossover regime has 
been intensively studied [l|, [2, |3, 0, d, @|. In an ultra- 
cold atomic system, the p-wave Feshbach resonance was 
also found, which allows one to tune p-wave interactions 
between atoms in the spin-polarized fermi system 0. It 
has now become possible to produce p-wave molecules 
between "^^K atoms using this technique Q. 

Atomic boson-fermion mixed systems, too, are ex- 
pected to show many interesting phenomena. In this 
mixed system, it has been reported the observation of 
simultaneous quantum degeneracy of Bose-Einstein con- 
densate (BEC) and Degenerate Fermi gas tS^iJAilli], and 
the collapse of fermions in the attractively interacting 
mixture |T^ . \\^ . On the theoretical side a number of 
studies on the static and d yna mical properties of the mix- 
ture have been made [IlIlE [H, [13, [ll] . We note here 
that the strongly coupled boson-fermion pair may behave 
as a heteronuclear molecule or resonance, and its role in 
the mixture has been studied [H, 113, [Hi, [H, [H. Re- 
cent experiments show an existence of the Feshbach res- 
onances between bosons and fermions [l^, [2^ , and quite 
recently a formation of the boson-fermion heteronuclear 
molecule in the optical lattice has been reported [1^ . As 
in the case of fermionic system, the existence of the Fes- 
hbach resonance allows one to control the boson-fermion 
interaction [27l . [28j . Tuning the interaction one may 
induce a collapse of fermions or a phase separation of 
bosons and fermions, depending on the sign of the inter- 
action [ii, [H, [13, [1|. 

In many experiments of atomic Bose-Fermi mixture, 
all atomic fermions have the same spin components as 



they are trapped by magneto-optical trap together with 
bosons. The fermion-fermion interaction has been neg- 
ligible since the s-wave scattering amplitude between 
fermions vanishes because of Pauli exclusion principle, 
and the other higher partial waves do not contribute at 
low temperatures. In a boson-fermion mixed system, bo- 
son density fluctuations may give rise to an attractive 
interaction between fermions as in the electron-phonon 
system, and may induce a superfluid transition [H, [s^l ■ 
This mechanism has been studied in the liquid ■^He-^He 
mixtures ^31.] . If this mechanism is strong enough in the 
ultracold atomic mixture, p-wave Coop er p airs will be 
formed between spin-polarized fermions [sj [H, [s^l • We 
should note, however, that a realization of the super- 
fluid depends on a balance of different mechanisms as 
mentioned above. A sufficiently strong attractive boson- 
fermion interaction, for instance, may cause a collapse 
of the system rather than the superfluid state. On the 
other hand, a strong boson-fermion repulsion may favor 
a phase separation of bosons and fermions instead of the 
superfluid transition. 

The purpose of the present paper is to investigate 
a possibility of p-wave superfluid transition in the re- 
pulsively interacting spin-polarized Bose-Fermi mixture 
at zero temperature. We consider superfluid transi- 
tion induced by a Bogoliubov phonon-mediated fermion- 
fermion interaction, and compare energies of the p-wave 
superfluid state and the boson-fermion phase-separated 
state. In view of the recent development in the tuning of 
boson-fermion interaction via Feshbach resonances, we 
study the system in a broad range of the interaction 
parameters. Earlier study in this direction shows that 
the p-wave superfluid transition is hard to occur, i.e., its 
transition temperature is too low to be attainable, be- 
cause the phonon-mediated attractive interaction is very 
weak [32| . A stronger repulsive interaction would cause 
an instability of the system towards phase separation. 
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and then the induced attractive interaction reduces. We 
show, however, that for a very strong repulsive interac- 
tion it is possible to realize the p-wave superfluid tran- 
sition. In the superfluid phase, anisotropic energy gap 
appears like that of p-wave superfluid of '^He and heavy- 
fermion systems. We also show that the Cooper pair is a 
tightly bound state like a diatomic molecule in the strong 
boson-fermion coupling regime 35, 36]. 

Main content of the present paper is as follows. Sec. II 
derives ground state energies of the system in the p-wave 
superfluid phase due to the phonon-mediated attractive 
interaction and in the normal phase by the method of 
Ref.^29i]. Sec. Ill presents phase diagram of the system 
by direct energy comparison of the p-wave superfluid and 
phase-separated states. The momentum dependence of 
the energy gap and other properties of the superfluid 
state are shown. Finally, Sec. IV is a summary. 



experimental situation [U [2^, [23, . Hereafter, we take 
obf > 0, leading to repulsively interacting Bose- Fermi 
mixtures. 

We consider two types of quantum phases of fermions, 
that is, the normal phase and the superfluid phase. On 
the one hand, a phase separation of bosons and fermions 
occurs in the normal state of fermions, as discussed in 
Ref. [1^ for a large positive Ubf- On the other hand, 
for a strongly repulsive boson-fermion coupling, an effec- 
tive fermion-fermion interaction via density fluctuations 
of BEG becomes strong and is expected to lead to a su- 
perfluid state of fermions. In the following we calculate 
ground state energies of these two types of phases, sepa- 
rately [33. 



A. Superfluid phase : p-wave pairing 



II. MODEL 

We consider a spin-polarized uniform mixture of 
atomic bosons of mass me and fermions of mass mp at 
zero temperature. The system is described by the Hamil- 
tonian 



H — Hp + Hb + Hi, 



where 

Hp = X!('^k-MF)4ck: 

k 



(1) 



(2) 



E^l+q^k-q^fck, (3) 



k kpq 

Htnt = Ubf ^ &p+qcLq^pCk- 
kpq 



(4) 



Here c^"^ and 6^^^ are the annihilation (creation) opera- 
tors for the fermionic and bosonic atoms of momentum k, 
respectively. The corresponding momentum states have 
kinetic energies of bosons = h^k^ /2mB and of fermions 
ejj^ — h^k^ /2mF. The chemical potentials of fermions and 
bosons are denoted by and fiB ■ The boson-boson and 
boson-fermion collisions are described by the interaction 
strengths Ubb = ^Trfi^ ubb / mB and Ubf = ^.-kH^ ubf I riir, 
respectively, where obb and obf are the corresponding 
s-wave scattering lengths and = mBm-p/{mB + m-p) 
is the reduced mass. The elastic fermion-fermion s-wave 
scattering for spin-polarized fermions is absent because 
of Pauli exclusion principle. 

In order to give a realistic estimate of the phase di- 
agram, we consider a system of ^^Rb-'^^K mixture for 
the atomic bosons and fermions, where me = 1.419 x 
10-25 kg and mp = 0.649 x lO'^^ kg, and obb = 98.98ao 
with Bohr radius cq. As for the boson-fermion interac- 
tion, we assume the scattering length to be tunable via 
Feshbach resonance. This can be realized in the current 



We assume that all bosons are condensate and the fluc- 
tuations are treated by excitations of Bogoliubov phonon 

of energy: huj^ ~ e^(Sk + 2C/bb"-b) where tib is a con- 
densate density. Thus the bosonic part of the Hamilto- 
nian, Eq.(3), is transformed to 

Hb = -MB^B + C/BB^^|/2 + ^ riioiplp^, 

k#0 

and the boson-fermion interaction Hamiltonian, Eq.(4), 
is also transformed to 



Hint — UBFnB E c|^Ck 



%F\/n^E 



k#0,p 



k+p'^p 



Pi 



k'^p'^k4 



(5) 



where /J^^-* is a Bogoliubov phonon annihilation (creation) 
operator. In these transformations, we have neglected 
higher order terms of (3'^\ We regard the flrst term of 
Eq.® as an energy shift of the fermionic chemical po- 
tential /ip = /ip — C/bf^ib- 

By applying the second-order perturbation theory, the 
phonon-mediated interaction between two fermions is de- 
rived as one-phonon exchange process. 



Hi, 



E 

pp'q 



t^FBF(p',q)cj3+q( 



with 



C/FBF(p',q) = Uby^B 



2e^ 



We assume that Fermi velocity up = hk^/m-p is 
much smaller than Bogoliubov sound velocity s = 



3 



\/ U-q-qu-q I m-Q , i.e. vy <C s, so the phonon-mediated 
effective interaction can be written as in Ref . 13211 . 



UvBFiq) 



1 



C/bb 1 + {h\q\/2mBs) 



(6) 



This interaction is equivalent to an effective interac- 
tion between fermions induced by density fluctuations 
of background of BEC when the retardation effect is ne- 
glected [29, 30]. The effective Hamiltonian of fermions is 
then given by 



-\Y. C^FBF(k-k') 



^''p/2+k'^P/2-k'^P/2-k'Cp/2+k'- 



Pkk' 



(7) 



Since we consider spin-polarized fermions, the effective 
interaction in the channel with even angular momentum 
I is absent due to antisymmetrization of the orbital wave 
function in the relative coordinate. Furthermore for the 
effective interaction described by Eq. ^ the contribution 
to the interaction for higher I can be negligible [30| • Thus 
we extract the dominant / = 1 component of the phonon- 
mediated interaction 



(8) 



i—x,y,z 



with 



+ {2mBs/hY 



Akk' 



In 



{k + k'Y + {2mYis/hY 



{k-k'f + (2mBs/n)2 



where k = |k| and ki = (k/|k|)i. In this approximation, 
only p-wave pairing is possible to realize. 

Let us consider a Cooper pair with zero center-of-mass 
momentum and introduce a p-wave pair energy gap as 

A(k) = E[/|;Bp(k,k')(c_k'Ck') 



where ( ) denotes an expectation value in the ground 
state. In the standard BCS theory the effective Hamilto- 
nian Eq. ([7]) can be diagonalized by Bogoliubov transfor- 
mation Q^k — UkCjj. — WkC^k' 



k 



|A(k) 
2-Bk 



, (9) 



where the quasi-particle energy is E'k = ■\/^^ + |A(k)P 
with ^k = ^k ^ /^Fi ^^"i ""k and Uk are given by 
Mk = V(l + Ck/£^k)/2 and Dk = v/(l - Ck/£;k)/2, re- 



spectively. The gap equation 

A(k) = ^{/|;3p(k,k')^, 



(10) 



and the requirement of the mean number of fermions 



n-p 



determine the energy gap A(k) and chemical potential of 
fermions /ip . Concrete forms of the gap will be discussed 
in the next section. 

Once we have the optimized value of the energy gap 
and chemical potential, we obtain the ground state en- 
ergy per volume in the superfluid phase as 



2 ^ V 2£;k 

/^f'^F + ^C/bB"-b + %F"-B?lF- (12) 



ESF 



The third and fourth terms in Eq. p2)) are mean-field 
energies of the boson-boson and boson-fermion interac- 
tions. 



B. Normal phase : phase separation 

Now we consider a normal phase of repulsively interact- 
ing Bose- Fermi mixtures. In Ref. (29| . Viverit et al. stud- 
ied a phase diagram of the uniform mixtures at zero tem- 
perature and predicted three types of equilibrium states: 
(A) a single uniform mixed phase, (B) a purely fermionic 
phase coexisting with a mixed phase, and (C) a purely 
fermionic phase coexisting with a purely bosonic one. 

By the method to find equilibrium states in Ref. [29l |. 
we obtain the ground state energy of the system of a finite 
volume V and total number of bosons A^b and of fermions 
Ny- Suppose that the system is composed of two phases 
of volume Vi and V2 where V = Vi + V2. The numbers of 
bosons and fermions are given by N^i and N-pi for z-th 
phase, respectively, yielding to Nb — Nbi + Nb2 and 
Np = Nfi+Nf2- Thus the total densities of bosons ub — 
N-q/V and of fermions np = Ny/V are given by tib = 
nBi'y+nB2(l— '!') andnp = nFiw+nF2(l— w), respectively, 
where v = Vi/V. Total energy per volume V is given by 



En = Eiv + E2{l-v), 
where the energy per volume Vi of i-th phase has 



(13) 



3 1 

Et ^ -f-FiUFi + -UBSn^i + [/BF"Bj"-Fi- (14) 

5 2 



In Eq. (|M)) . the first term is the kinetic energy of 
fermions where epi = h'^{6TT'^nFiY^^ /2mF is Fermi en- 
ergy of i-th phase. The second and third terms cor- 
respond to the mean-field energies of the boson-boson 
and boson-fermion interactions. In j-th phase, the pres- 
sure is Pi = —d{EiVi) / dVi and the chemical potentials 
of bosons and fermions are given by /iBi — dEi/dn^i 
and ^Fi — ^Fi — dEi/duFi, respectively. The equilib- 
rium conditions of the pressure and chemical potentials 
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of bosons and fermions between two phases, Pi = P2, 
A*Bi = ^62) and epi = €^2, determine Vi, riBi, and npi- 
As a result, the three types of equihbrium states can be 
reahzed. 

Figure [ija) shows a phase diagram for the ^^Rb-'^^K 
mixed system in the hf/ub vs. Ubf/Ubb space for 
ns = 10i4cm-3 and Ubb = 5.157 x IQ-^^ [Jm^]. The 
lower and upper solid lines in Fig. [Ija) represent the 
boundaries between the phases (A)-(B) and phases (B)- 
(C), respectively. The result shows that a phase sepa- 
ration of bosons and fermions is preferred for a higher 
fraction of fermions and for a stronger boson-fermion 
interaction strength compared to the boson-boson one. 
The dashed line corresponds to the critical curve above 
which the uniform mixture is unstable against small den- 
sity fluctuations. 

Figure [ijb) shows total energy per volume of the 
ground state En calculated from Eqs. ([T^ and as a 
function of Ubf/Ubb for np — ub — 10^^ cm^'^. The to- 
tal energy per volume of the ground state scaled by epup 
with ep = ?i^(67r^nF)^/^/2mF is plotted against the inter- 
action ratio i7BF/t^BB- The energy value increases as the 
repulsive boson-fermion interaction becomes stronger un- 
til Ubf/Ubb — 6 above which fermions start to separate 
from the bosonic cloud. For Ubf/Ubb }t 8 the ground 
state is the completely phase-separated phase (C), and 
the energy becomes independent of the interaction ratio 
Ubf/Ubb- 



III. RESULTS 

In this section, we examine phase transition between 
normal state and superfluid state of fermions. To do 
so, we assume two typical types of energy gaps, an ax- 
ial state and a polar state in the superfluid phase, and 
calculate ground state energies in the two phases and 
directly compare them. We obtain phase diagram of 
^^Rb-^°K mixtures showing that phase transition from 
phase-separated state in the normal phase to p-wave su- 
perfluid state occurs as repulsive boson-fermion coupling 
becomes stronger. The unique characters of the p-wave 
pairing, such as momentum dependence of energy gap, 
momentum distribution of fermions, and quasi-particle 
excitation spectrum, in a strong boson-fermion coupling 
regime are also discussed. 



A. Transition from phase separation to p-wave 
superfluidity 

In order to obtain ground state energy of superfluid 
fermions, we need to calculate the total energy described 
by Eq. (|12p . In general, the p-wave pair energy gap can 
be expanded by spherical harmonics 
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A(k) = Yl A'"(A:)yi,™(0, 
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Pure boson 

and pure fermion phase 
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FIG. 1: (a) Phase diagram of normal phase of the ^^Rb-^^K 
mixtures. There exist three types of equilibrium states: (A) 
a single uniform mixed phase, (B) a purely fermionic phase co- 
existing with a mixed phase, and (C) a purely fermionic phase 
coexisting with a purely bosonic one. The lower and upper 
solid lines in Fig. 1(a) represent the boundaries between the 
phases (A)-(B) and phases (B)-(C), respectively. The dashed 
line corresponds to a critical curve above which the uniform 
mixture is unstable against small density fluctuations, (b) 
Total energy per volume of the ground state En scaled by 
epTiF as a function of Ubf/Ubb for np = ns = 10^*cm~^. 



We assume, however, just two typical types of pair energy 
gaps, so-called an axial state [3g| 



A-'^(k) = A-'^(fc) it, + iky), 



and a polar state 



AP'(k) = AP^(fc) k 



(15) 



(16) 



and regard a lower energy state of them as the ground 
state in superfluid phase. Note that there is no symmet- 
rical solution due to the spin-polarization of fermions in 
this system. This fact is in stark contrast to unpolar- 
ized Fermi systems such as '^He and heavy- fermion sys- 
tems [11] . Amplitudes of the axial and polar states have 
|A'^'=(k)| = A^'^(fc) I sin 6*1 and |APi(k)| = AP\k) \ cose\, 
corresponding to anisotropic energy gaps in momentum 
space with zero energy gap on the north and south poles 
for the axial state and on the equator for the polar state, 
respectively. 

Figure[2)Ja) shows total energies per volume of the axial 
state (tetragon) and polar state (asterisk) in unit of eprip 
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FIG. 2: (a) Total energy per volume of the system Esf 
scaled by eprip and (b) chemical potential fip scaled by ep as 
a function of Ubf/Ubb at np — ub = lO^^cm"^. Each mark 
of a tetragon and an asterisk denotes the axial state and polar 
one. 



as a function of Ubf/Ubb for np — ub — 10 cm"^. It 
shows that the total energies per volume of both of axial 
and polar states increase as the boson-fermion coupling 
becomes stronger until Ubf/Ubb — 30. It is caused by 
the mean-field energy of boson-fermion repulsive inter- 
action. At sufficiently strong boson-fermion interaction 
strength, however, the effect of the chemical potential 
is significant and the total energy turns to decrease. We 
plot the behavior of the chemical potential of the super- 
fluid states function of Ubf/Ubb in Figure[2l[b). 
In the weak boson-fermion coupling regime, the relation 
/ip ~ ep holds as in the ordinary BCS theory. However 
in the strong boson-fermion coupling regime, fip starts 
to decrease, eventually becoming negative with increas- 
ing Ubf/Ubb- Negative values of the chemical potential 
indicate formation of bound pairs [ssl . [36j . 

The dot-dashed line in Fig. [2{a) corresponds to total 
energy per volume of the ground state in normal phase, 
which is the same as the result in Fig. [IJb) . The direct 
comparison between the ground states in the two phases 
shows that phase transition from the phase-separated 
state in normal phase to the axial state in superfluid 
phase occurs when the interaction ratio Ubf/Ubb in- 



FIG. 3: Phase diagram with considering the superfluid transi- 
tion at riB = lO^'^cm"'^. Circles denote the boundary between 
the p-wave superfluid phase with axial state and normal phase 
with phase separated state. 



creases. We see that the energy of the axial state is lower 
than that of the polar state in the strong boson-fermion 
coupling regime. It is caused by the angle dependence 
of the energy gap. The axial state has anisotropic en- 
ergy gap in momentum space with zero energy gap on 
the north and south poles. For the polar state, on the 
other hand, the energy gap vanishes on the equator in 
momentum space. So the former state contributes to 
the decrease of the total energy more than the latter 
case. Figure [3] shows phase diagram of the system in the 
Ubf / UBB-nF / nB plane for ub = lO^^'cm"'^. The lower 
and upper lines are boundaries of the phases (A)-(B) 
and the phases (B)-(C) in normal phase as already ex- 
plained in Fig.[lja). The circle corresponds to the critical 
points below and above which the phase-separated state 
and axial state realizes, respectively. The dependence on 
np/uB of the phase boundary is not strong. Thus we 
conclude that p-wave superfluidity can be observed at 
sufflciently strong boson-fermion interaction compared to 
boson-boson interaction. This is the main result in the 
present paper. Using typical experimental parameters, 
we estimate the s-wave scattering length between boson 
and fermion simply from the rate of coupling constant. 
It is estimated as aBp/aBB — 23 that is tunable by the 
technique of Feshbach resonance [13, ■ 



B. p-wave pairing state in strong boson-fermion 
coupling regime 

We have just seen that p-wave superfluid state arises 
at sufflciently strong boson-fermion interaction strength. 
Properties of the strong boson-fermion coupling super- 
fluid state such as energy gap, quasi-particle, and mo- 
mentum distribution are expected to differ from those of 
the weak boson-fermion coupling one. In the following 
we will show some aspects of the p-wave pairing in the 
strong boson-fermion coupling regime. 

Figure [5] shows the radial component of energy gap 
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FIG. 4: Momentum |k| dependence of the radial gap of the 
axial state A^''(fc) scaled by ep at np = ne = 10^*cm~^. 
Solid, dashed, dashed-dotted, and dotted lines denote results 
for t/Bp/C/BB = 20, 30, 40 and 50, respectively. 



FIG. 6: Population numbers of fermions, v^. — (1 — ^k/£^k)/2, 
in momentum space. Solid, dashed, dashed-dotted, and dot- 
ted lines denote Ubf /Ubb = 10, 20, 40 and 50, respectively. 




FIG. 5: Quasi-particle energy in the axial state _Bk scaled by 
ep as a function of fc/fcp with the angular component — n/2 
at np = nB = lO^'^cm"^. Dashed, dashed-dotted, and dotted 
lines denote ?7Bp/t^BB = 20, 40, and 50. And a solid line 
denotes dispersion relation of ideal fermi gases. 



of the axial state A'^^{k) in unit of Fermi energy ep as 
a function of k/kp where kp — (GTr^np)^^'^. The energy 
gap peaks at around the Fermi surface, and a Cooper pair 
can be formed even in high momentum region of |k| > fcp. 
The gap also becomes larger, and the width around the 
peak gets broader with increasing Ubf /Ubb- The gap 
vanishes at A; = 0, and it originates from the property of 
the interaction with angular momentum I = 1 between 
fermions. 

Since A(k) corresponds to a pairing potential that af- 
fects on the pairing state of (k, — k), the gap crucially 
influences quasi-particle spectrum, see Figure [H The 
quasi-particle energy scaled by Fermi energy ep is plotted 
as a function of fc/fcp. The solid line denotes dispersion 
relation of a free fermion. In the weak boson-fermion cou- 
pling regime, the quasi-particle energy differs from that 
of a free fermion only around the Fermi surface. In the 
strong boson-fermion coupling regime, the quasi-particle 
spectrum changes dramatically. For /ip < 0, the gap 
opens at zero momentum and corresponds to molecular 
binding energy [H, . We can see that a Cooper pair 
is now a tightly bound state of Z = 1 like a diatomic 



molecule with the binding energy — 2|/ip| in the strong 
boson-fermion coupling regime. 

Finally, the momentum distribution function of 
fermions for different Ubf /Ubb values is shown in Fig- 
ured As Ubf /Ubb increases, the occupation of fermions 
reduces significantly and Fermi statistics turns to be 
less important. The momentum distribution at large 
Ubf /Ubb is expected to be proportional to momentum 
distribution of ^ = 1 bound state of two fermions. The 
behavior of the distribution at fc = is explained as fol- 
lows. From Eg. ljlip the population number of fermions 
at A: = becomes ^^=0 = (1 + Mf/IMfD/S- In the 
Ubf /Ubb = 10, 20, and 40 cases, the population num- 
ber of fermions at fc = becomes ^^=0 ~ ^ because of 
/^p > shown in Fig. [^l^b). In the Ubf/Ubb — 50 case 
with /Zp < 0, however, the population number of fermions 
vanishes at fc = 0. This result can be also understood by 
noting the following fact. For /j,p > 0, as indicated by 
Fig. m fermions of extremely small momentum are al- 
most non- interacting and occupy as free particles. For 
/ip < 0, the population number distribution tends to 
have the form of momentum distribution of single bound 
pair wave function [36l | and vanishes at fc = because of 
I = 1 bound state. The particle distribution is observ- 
able in the real experiments, for example, by using the 
time-of-flight method. Thus we consider this result as an 
indication for the observation of the p-wave superfluid. 



IV. SUMMARY 

We considered a system of repulsively interacting Bose- 
Fermi mixtures of spin polarized uniform atomic gases at 
zero temperature. We investigated the possibility of real- 
ization of p-wave superfluidity of fermions due to effective 
attractive interaction via density fluctuations of BEC in 
the case of strong boson-fermion interaction within mean- 
field approximation. By direct energy comparison be- 
tween p-wave superfluid and phase-separated states, we 
found that p-wave superfluidity can be observed at suf- 
ficiently strong boson-fermion interaction compared to 
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boson-boson interaction. 

We also discussed unique features of the ground state 
in strong coupling p-wave superfluidity. We calculated 
the quasi-particle energy and found that a Cooper pair 
is a tightly bound state like a diatomic molecule in the 
strong boson-fermion coupling regime. We also calcu- 
lated the momentum distribution function of fermions. 
In the strong boson-fermion coupling regime, we showed 
the property of the single bound pair wave function and 
suggested an observable indication of p-wave superfluid 



in the real experiment. 

The present analysis was made by using the standard 
mean- field approach with BCS wave function [ssf. The 
results at very strong coupling would be modified by 
fermionic self-energy correction as well as by dynamical 
screening of the effective interaction |4(| . The consider- 
ation of these effects is left for a future study. Neverthe- 
less, the mean-field approaches are expected to provide 
good ways of building up intuition. 
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